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Abstract 

The microscopic origin of dissipation of a driven quantum many body system 
is addressed in the framework of a parametric banded random matrix approach. 
We find noticeable violations of the fluctuation-dissipation theorem and we observe 
also that the energy diffusion has a markedly non-Gaussian character. Within the 
Feynman- Vernon path integral formalism and in the Markovian limit, we further 
consider the time evolution of a slow subsystem coupled to such a "bath" of intrinsic 
degrees of freedom. We show that dissipation leads to qualitative modifications of 
the time evolution of the density matrix of the slow subsystem. In either the spatial, 
momentum or energy representation the density distribution acquires very long tails 
and tunneling is greatly enhanced. 
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1 Introduction 



It is by now an accepted fact that spectral fluctuations of many quantum systems are 
identical to those resulting from random matrix theory. Random Hamiltonians have 
initially been introduced by Wigner to describe nuclear properties at excitation energies 
of the order of the nucleon binding energy (T| . They have found applications in many areas 
of physics. Just to cite several directions, they have been used to describe properties of 
simple quantum chaotic systems @, spectra of complex molecules ||, condensed matter 
physics 0, nuclear properties and spectral fluctuations obtained in QCD calculations 
||. There is enough experimental evidence from all these physical systems to justify the 
assignment of the qualifier universal to their spectral properties. 

In this article we address the problem of the evolution of a quantum mechanical system 
for which the spectral fluctuations are assumed to be universal, the average density of 
states has the proper thermodynamic behaviour and which, in addition, is capable of 
changing its "shape". By "shape" we understand either some parameters describing the 
actual geometrical shape of the system or some other similar global characteristics, or 
externally applied fields. The shape can either be controlled by the experimenter, or can 
have its own quantum dynamics. Some classic examples of the first type of system are 
those studied in thermodynamics, while examples of the second type are atomic nuclei 
and complex molecules. We are interested in what happens when the shape changes at a 
finite rate. In other words, we are interested in the microscopic nature and character of 
dissipation in finite quantum many-body systems. 

The paper is organized as follows. In Section II we describe the structure of the 
Hamiltonian and the evolution equations for driven systems and what we often refer to 
as a simple quantum system in interaction with a complex environment. The theoretical 
framework in this last case was developed many years ago by Feynman and Vernon |J 
and used extensively, particularly in condensed matter physics, in Caldeira-Leggett type 
of treatments of dissipative tunneling Q. There also exists in the literature another type 
of approach, somewhat more phenomenological in nature, namely the quantum Langevin 
equation ||, which is claimed to lead to a dissipative dynamics similar to that of the more 
popular Caldeira-Leggett formalism. Our main contribution consists in the introduction 
of a new type of environment or "bath", based on a parametric banded random matrix 
approach ||, IU|, |TI |. In our opinion, such an approach, besides being extremely flexible, 
has a better microscopic foundation and is better suited to describe finite many-body 
quantum systems. 

In Section III we study the temporal evolution of driven complex quantum systems 
which we describe using a parametric banded random matrix approach. The shape is an 
externally modulated parameter. Should this shape change infinitely slowly, the evolution 
of the system would be reversible. However, at any finite rate of change, it is for all 
practical purposes irreversible. A simple example should suffice. Assume that we have 
a pump, inside of which the motion is chaotic. In the pump there is one single particle, 
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bouncing elastically off the wall and off the piston. At some point in time someone starts 
moving the piston, eventually bringing it back to its initial position. We assume that 
during the entire time the piston is in motion, the person is not permitted/allowed or 
able to get any information about the position or momentum of the particle. Further, 
the person is not allowed to act in any direct way on the particle. For almost all closed 
trajectories of the piston, covered with a finite velocity, the particle inside will not return 
to its initial state. One can then classify the "state transformation" of the particle in the 
pump as irreversible. This in some sense paraphrases the almost century old argument 
between Boltzmann on the one side and Zermelo and Loschmidt on the other In 
that argument the role of the particle inside the pump was played by the infinite number 
of atoms in a gas. Part of the Boltzmann's argument was that because their number 
is infinite, one cannot in any conceivable way actually reverse their velocities, and thus 
irreversibility arises. We have only changed "cannot" with "is not permitted" and arrived 
at the same result, but now for perhaps one of the most simple systems possible. This is 
essentially what one often encounters in real situations. Even with the intrinsic system 
having a finite number of degrees of freedom, there is no obvious way by which one can 
control their state or acquire information about their microscopic state without disturbing 
them, and in particular reverse their momenta. One can only control some "external" 
parameters. 

In Section IV we analyse a few cases of a simple quantum system, coupled to such a 
complex, but finite, "environment" within the adiabatic approximation. In this case the 
energy transfer is allowed in both directions, from one to the other subsystem and back as 
well. However, if one of the subsystems has a large number of active degrees of freedom 
the energy flow will most of the time occur in one direction. This is a generic situation, 
encountered in numerous quantum finite many-body systems. 

A short summary of our results and an outlook for future investigations are presented 
in the last Section V. 



2 Evolution equations 

The Hamiltonian governing the dynamics of a quantum system coupled to a complex 
environment is assumed to be of the form 

H(X,x) = H (X)+H 1 (X,x). (1) 

We refer to X as shape variables. In the ensuing formulas we shall not display explicitly the 
dependence of the Hamiltonian on the intrinsic variables x of the environment, but rather 
discuss its matrix elements in a fixed intrinsic basis. The part of the total Hamiltonian 
(H) which depends on the intrinsic coordinates H\(X) is defined as a parametric banded 
random matrix, whose matrix elements depend on the "slow" coordinate X 

[H 1 (X)] fj = [h ) ij + [h 1 (X)] ii . (2) 
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ho is taken to be diagonal and defines the average density of states, with (k\ho\l) = 
[ho]ki = £k8ki- We refer to these eigenstates as "typical states" of the intrinsic system 
with an energy e. In Refs. ||, [R], |TT| we have discussed at length the reasons why one 
chooses this specific form of the Hamiltonian. For an intrinsic subsystem with a large 
number of degrees of freedom the average density of states, 



P (e) = Ti5(H 1 (X)-e) 



(3) 



for each given shape X increases sharply with energy. The overline denotes here a proce- 
dure for extracting the smooth part of p(e) as a function of energy and it amounts essen- 
tially to an ensemble average, to be introduced below. For a many Fermion system, p(e) 
has a roughly exponential behaviour. Recall that In p(e) is approximately proportional to 
the thermodynamic entropy of the intrinsic system, which is an extensive quantity. The 
fact that the average density of states for the intrinsic subsystem has such a behavior 
is a key element of the entire approach. This is equivalent to stating that the intrinsic 
subsystem has a large heat capacity and thus can play the role of a "reservoir" , although 
not necessarily ideal. In principle p(e) can be X-dependent as well, but we shall ignore 
this aspect here. Without an X-dependence of the average density of states, mechanical 
work cannot be performed on or by the model environment we study here, and only heat 
exchange is allowed. 

In the basis of the eigenstates of h , we define hi(X) as a parameter dependent N x N 
real Gaussian random matrix, which is completely specified by its first two moments 



Mx)} kl 



o. 



[M*)k#iQ0]ju = Wji + SMQijiX-Y). 



(4) 



The overline stands for the statistical average over the ensemble of random Gaussian 
matrices. Even though we shall limit all the formulas in this work to the Gaussian 
orthogonal ensemble (GOE), with minor changes, the formalism is equally applicable to 
other Gaussian ensembles. In Eq. ([|), Qij(X — Y) is a "bell shaped" correlation function 
with a characteristic width Xq. Physically this is simply the formal expression of the 
fact that intrinsic states corresponding to vastly different values of the shape X have 
statistically independent compositions. The dependence on i,j allows for the description 
of banded matrices, where an effective number of states N < N are coupled by hi(X). 
It is convenient to use an explicit parameterization of Qij, which explicitly incorporates 
the average density of states and the bandwidth of the statistical fluctuations[|l3|]: 
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(5) 



Here G(x) = G{—x) = G*(x) < 1, G(0) = 1 and is the spreading width for the 
intrinsic subsystem. k (linked with the effective band width N Q w kqp(e)) and X are 
also characteristics of the intrinsic system. This parameterization is consistent with the 
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Figure 1: A portion of the instantaneous eigenvalue spectrum E n (X) as a function of the 
shape X. The Hamiltonian is defined through Eqs. (0), (|]) and (H), with unit average 
density of states and = 2rr, k = oo, G(x) = exp(— x 2 /2) and X = 1. 

experimental evidence that in many-body systems the spreading width changes rela- 
tively slowly with the excitation energy [|i~3|1 . Moreover, distributions of matrix elements 



extracted from various theoretical many-body models show strong deviations from a pure 
Gaussian distribution |T3| and are in qualitative agreement with Eq. (|5]). A typical Born- 
Oppenheimer spectrum (with an average unit level density) as a function of the shape X 
is shown in Fig. 1 for a particular realization of the random Hamiltonian H\(X). For 
each fixed value of the shape variable X the spectrum is characterized by fluctuations 
very similar to GOE. 

2.1 Driven subsystem 

A particular situation of definite physical interest is that when the time dependence of 
the "slow " variables is known and/or controlled externally as for example in the case 
of an applied external electric or/and magnetic field(s). This is a typical situation in 
thermodynamics, when the controlled parameters are changed adiabatically, insuring the 
reversibility of the transformation. At any finite rate of change the transformation looses 
its reversible character and dissipation sets in. In an experiment, if one has some control 
over the parameters, as a rule, one has little control on the intrinsic subsystem. Thus, if 
we move the parameter X(t) at finite velocity along a closed circuit in parameter space, 
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returning to the initial point, the final intrinsic states will in general be distinct from the 
initial one. These changes of the intrinsic subsystem, which not only depend on the rate 
of change of the controlled parameters, but also on the specific paths in the parameter 
space, are referred to as dissipative effects. State changes of the intrinsic subsystem, 
which depend on the path in the parameter space, but do not depend on either the 
rate at which this path is covered or on the direction in which the path is travelled, 
can be incorporated into reversible type of transformations, by slightly generalizing the 
definition of adiabatic transformations Such effects can be linked with the appearance 



of effective abelian or nonabelian gauge fields besides the familiar (in the Born- 
Oppenheimer approximation) effective potential forces. 

For the sake of simplicity we assume that X(t) = V^t and that the average level 
density has an exponential form: p(e) = poexp((3e). Thus (3 — 1/T — d In p(e)/de can be 
interpreted as the inverse thermodynamic temperature of the intrinsic subsystem. (Note 
that the spectrum extends to infinity in both directions from e = 0.) 

The time evolution of the fast subsystem is found by solving the time-dependent 
Schrodinger equation in the form [[], |TTfl : 



4>(t) = T exp 



-- fdsH^s)) 
n Jo 



0(O)=W(X(t))0(O). 



(6) 



where T is the time ordering operator, and U(X(t)) the propagator. (We assume that the 
initial state 0(0) is uncorrelated with the Hamiltonian Hi(X(t)) at later times; correlated 
initial conditions have been discussed elsewhere One can show that in the leading 
order in an expansion in 1/N the average propagator U(X(t)) = U(X(t)) is diagonal in 
the representation we have chosen. Its diagonal matrix elements have the following form 



U k {X(t)) = (A;|Texp 



-ijf dsHi(X{8)) 



\k) = ^[--f-)a(X(t)) (7) 



( note that a(X(t)) is state independent) and a(X(t)) satisfies the following integral 
equation |fLl] : 

si — s 2 )V \ , s 



a(X(t)) = 1 - y Jd Sl J% 2 a(X( Sl - s 2 ))a(X(s 2 ))P( Sl - s 2 )G f- 
Here P(s) is given by 



P(s) = P*(-s) 



2nh 



exp 



s + % 



(9) 



when the correlator [hi(X)]ij[hi(Y)] k i is defined as in Eq. ([5]). In order to be able to 
compute averages of observables, we need to introduce the set of generalized occupation 
number probabilities 



AT fc (X(ti),X(t 2 )) = (0(ti)|*)<A#(t a )) = E {l\UKX{ti))\k){k\U{X{t 2 ))\l)m{Q). (10) 

i 
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Thus rikit) = J\fk(X (t) , X (t)) is the occupation probability of the state \k). The general- 
ized occupation number probabilities can be determined by solving an evolution equation 
for the characteristic functional 



AT(*(*i),*(f2),T) = (0(ti)|exp 



ih (r -ti + t 2 

h 



de k p{e k )M k {X{ti), X(t 2 )) exp 



ieJr-ti + t 2 ) 



The equal time functional J\f(X(t),X(t), r) is thus the Fourier transform of the occupation 
number probabilities n^t). J\f(t,t,r) is an extremely useful quantity, since it provides 
the following cumulant expansion ||19||: 



Af(X(t),X(t),r) = exp 

k 



V h 



exp 



„ n n! 



(12) 



where ((0(t)|/iQ are the cumulants. If we assume that the initial occupation number 

probabilities are n (0) = 1 and nj(0) = n_; (0) = for / ^ (remember that the spectrum 
of Hi(X) is infinite in both directions), then we find that J\f(X(ti), X(t 2 ), r) satisfies the 



evolution equation [ITT 



jV(X(ti),X(t 2 ),r) = a*(X(t 1 ))a(X(t 2 )) + ^ f tl ds 1 [ t2 ds 2 M(X(s 1 ),X(s 2 ),r) (13) 

a jo jo 

x (si - s 2 )Vo X 



xP(si - s 2 - r)G 



a\X(t 1 -s 1 ))a(X(t 2 -s 2 )). 



2.2 The path integral approach 

When the shape variables X become dynamical variables, there is energy exchange be- 
tween the two subsystems, and the dynamics becomes more complicated. A formalism 
to tackle the case when both subsystems have to be treated quantum mechanically has 
been put forward by Feynman and Vernon ||. One can write the following double path 
integral representation for the density matrix of the entire system 

H(X,x,Y,y,t) = [ dX dY ^(X )r(Y ) f X(t) ~ X VX(t) [ ¥{t) ~ Y VY(t) (14) 
J Jx(o)=x Jy(0)=Y 

x exp{~[S (X(t))-S Q (Y(t))}} 
x (x|Texp 

where T and T a represent the time ordering and time anti-ordering operators respec- 
tively and So(X(t)) is the classical action corresponding to the Hamiltonian Hq(X). The 
particular form for the initial state wave function we have used here, namely 

V{X,x)=ij){X)4>(x) (15) 



% Jo 



dt'H^Xit')) 



T a exp 



h Jo 



df'H^Yit")) 



\y), 
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is not unique and other choices are equally possible, for example, a density matrix. By 
introducing the influence functional 



C(X(t),Y(t),t) 



T a exp 



ijT (y (*"))] }{Texp [-~jf dt'HiiXit')) 



(16) 

one readily obtains the following double path integral representation for the density matrix 
of the "slow" subsystem 



p(X,Y,t) 



f X(t)=X fY(t)=Y 

dX dY iP(X )i>*(Y ) / VX{t) / VY(t) 

JX(0)=X Jy(0)=Yo 



x exp \ - [S (X(t)) - S (X(t))] \ C(X(t), Y(t),t). 



(17) 



The formulation of the problem through a path integral representation serves only as a 
very convenient vehicle to obtain an evolution equation for the density matrix p(X,Y,t). 
(Our usage of the same Greek letter p for two different quantities, average level density 
for the intrinsic subsystem and density matrix for the "slow" subsystem, should not lead 
to confusion.) 



Using the formalism developed in Ref. |TTJ one can derive relatively simple analytical 
expressions for the influence functional. For the case of an adiabatic evolution of the slow 
subsystem, it was shown in Ref. |TTJ that the influence functional has the simple form 



C(X(t), Y(t),t) = M(X(t), Y(t), 0) = exp jf' [G(X(t') 7 Y(t')) - 1] dt' 



By combining the double path integral representation for the density matrix p(X, Y, t) 
with the above expression for the influence functional in the adiabatic approximation, one 
easily derives that the density matrix satisfies the following Schrodinger-like equation (for 
similar examples see Refs. 0]) 



ihdtp(X, Y, t) = {H {X) - H (Y) + T^ [G(X, Y) - l]}p(X, Y, t) 



(19) 



with the initial condition 



p(X,Y,0)=^(X)r(Y). (20) 
This equation for the density matrix describes a quantum mechanical Markovian process 



and it satisfies the conditions of the Lindblad's theorem 17 . Therefore the solutions 



of this evolution equation with meaningful physical initial conditions can be given a 
probabilistic interpretation at all subsequent times. In particular this means that 



p{X,X,t) > 0, 
dXp(X,X,t) = 1 



(21) 
(22) 



for any t > 0. 
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3 Temporal evolution of a driven system 



In this section we shall present analytical results for two limits, adiabatic and diabatic, 
and numerical results for the intermediate situation. In order to quantitatively determine 
whether the evolution of the driven subsystem (X(t) = VqI) is in either of these limits 
it is useful to introduce two time scales: i) the characteristic time scale for the slow 
motion r s i ow = X /Vo and ii) the characteristic time scale for the fast degrees of freedom, 
Tfast = H/hq. The adiabatic limit corresponds to t s i ow ^> Tf ast . (We shall also assume 
that in the adiabatic limit the condition k /3 <C 1 is also fulfilled.) The diabatic limit is 
obtained when T stow < r fast . 



3.1 Adiabatic limit 

The generalized occupation number probabilities J\f(X(t),X(t),r) can be computed by 
integrating Eqs. (§) and (0). The technical trick which is used is to replace in these 
equations the quantity P(s) introduced in Eq. @ with an appropriately chosen Dirac 



^-function |TT]. In the adiabatic limit P(t) is much narrower than G(Vot/X ). In the 



leading order in the parameter Kq(5 <C 1 one obtains for t > \r 

Wo 



M(X(t),X(t),r) = exp 




G 



X 



(*-M)--jP|. (23) 



In the strict adiabatic limit Vq — > 0, G(tVq/Xq) — > 1, and the first term in the exponential 
vanishes. The occupation numbers rik{t) reach rather quickly the asymptotic distribution 

1 r£ 

n k = o, (24) 

*4 + (W 



This Lorentzian shape is identical with the constant random matrix theory result || \L1 
During a time t w Tf ast , the slow variables hardly change and the dynamics of the fast 
system is almost identical to the dynamics governed by a constant random Hamiltonian. 
Our initial state in the middle of the spectrum, chosen as n (0) = 1, is thus spread over 
an energy interval ~ 1^ and the distribution has a Lorentzian shape. If the Hamiltonian 
is time independent, after this time there would be essentially no further evolution of the 
average occupation number probabilities. The subsequent dynamical evolution of the fast 
subsystem occurs only because the Hamiltonian Hi(X(t)) is time dependent, and only 
the subsequent time evolution of the system leads to dissipation and entropy production 
in the long time limit. 



One can now explicitly evaluate the cumulants ((</>(£) |/iQ|0(i))) |19| . All odd moments 



of ho vanish identically (since G(x) = G(—x) and thus there are only even powers of r in 
the expansion in Eq. ( p~2|) ). The reason for this is our assumption that k /3 — > 0, which 
we lift below. In the limit t — > oo, all even cumulants of ho increase linearly in time. If 
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G(x) = exp(— x 2 /2) (we shall use this form hereafter for illustrative purposes) then in the 
limit t — ■> oo 



D{V ) 



VH (W\ (2w)! 

(25) 









[2X 2 _ 




. ' slow . 



resulting in a non-Gaussian distribution. A Gaussian process would have only nonvan- 
ishing linear and quadratic cumulants. The energy diffusion constant is extracted from 
the time dependence of the second cumulant when t — > oo according to 



A|(t) = - E(t)}*m)) « (0(t)|^|0(t)) « const + 2£>(V )t, (26) 



where = (<f)(t)\Hi(X (t))\(f)(t)} . Note that the energy variance is time dependent only 
for a time dependent Hamiltonian. As the result of the symmetric initial distribution and 
since (3 — (remember that hq(3 <C 1), we obtain no friction (i.e. E(t) = const) and 
only nonvanishing even cumulants. To get friction we have to consider the next order 
corrections to the adiabatic limit k /3 <ti 1. Since for (3 > the average level density is 
increasing with energy, there will be on the average more transitions upward in energy 
than downward and thus the driven subsystem is heated up. One can show that the odd 



cumulants are then given by the following expressions [ 1 1 



((mhi n - l \m) = ^mmi n \m)- 

The case n — 1 corresponds to the Einstein fluctuation-dissipation theorem. In familar 
diffusive processes, all cumulants of order higher than the second are vanishing. The 
existence of large higher order cumulants results in energy tails of the energy distribu- 
tion significantly longer than in traditional phenomenological transport approaches, like 
Fokker-Planck or Langevin equations. One can show that in the tails the energy distri- 
bution has the following behaviour 

P(e) oc exp(— a\e\ hi 1/2 kl), (28) 

where a is some (time dependent) constant. A somewhat similar functional form has been 
determined for the distribution of conductance fluctuations in mesoscopic systems |2T| . 
The presence of these longer than expected tails is a clear indication that the excitation 
mechanism cannot be reduced to a simple random walk in energy space. 



3.2 Diabatic limit 

Another simple analytical solution can be obtained in the diabatic limit, when t s i ow = 
Xo/Vq Tjast = %/kq. In this case G(V s/X ) becomes much narrower than P(s) and 
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can be replaced with an appropriately chosen Dirac 5-function. The evolution equations 
(§) and ( pT3|) can be solved again and one obtains: 



Af(t, t, t) = exp 



K 



h 2 v 



K%(3^ 



exp 



-'Km: 



(29) 



(Note that in this case what we denote as the slow degrees of freedom are actually faster 
than the intrinsic ones.) This is similar to the functional form found in the adiabatic 
limit. Again, all the cumulants of ho increase linearly in time 



«0(t)|W(t)» 



h 2 V 6XP l 8 J lv/2, 



H. 



t. 



(30) 



where H n (a;) are Hermite polynomials, resulting in a non-Gaussian diffusion of the oc- 
cupation numbers. From Eq. (^) and the second cumulant, we find in this limit a 
completely different velocity dependence: 



D(V ) 



W ((3 2 k 2 + A) x / g/f 
8h 2 eXP \ 8 



1 



(31) 



The fluctuation-dissipation theorem, which provides the relation between the first and 
the second cumulants is in this case: 



(3D = 7 1 + 



[3 2 k 2 



(32) 



where jt = ((^>(t)\ho\^>(t))) . For a Gaussian diffusive process one would have obtained 
(3D = 7. Thus, in both adiabatic and diabatic limits we have obtained significant devia- 
tions from an excitation mechanism corresponding to a simple random walk in the energy 
space. One might expect that the intermediate regime will be very similar in this respect, 
as we shall substantiate in the next subsection. We thus conclude that the energy transfer 
to a "complex" quantum system is definitely not a Markovian process. The presence of 
the above mentioned long tails shows that the system retains for some time the memory 
of the direction it was proceeding and transitions into the same direction in energy are 
somewhat favored over transitions into the other direction. 



3.3 Arbitrary driving velocity 

For the intermediate velocity regime we have to resort to a numerical solution of the 
evolution equations Eqs. (W) and (|T3"D. In Fig. 2 we show Ukit) as function of the level 
number k. The tails have an almost pure exponential behavior, in agreement with the 
theoretical expectations, see Eq. fl28|) . In Fig. 3 we show the behavior of the diffusion 
constant, D(Vq), from the adiabatic to the diabatic limit for some values of the parameter 
(3. In all cases, -D(Vq) evolves from quadratic (in the adiabatic limit) to an inverse velocity 
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Figure 2: The time dependence of the occupation probabilities rik(t) (note the vertical 
logarithmic scale) plotted as functions of the level number k, for V = 4, (3 = and 
kq ps 20. The values of the remaining parameters are specified in the text. The narrower 
distribution corresponds to t — 1 with widening curves for t — 2, 3 and 3.5 respectively. 
The small notch in the uppermost curve is a remnant of the initial conditions n (0) = 1 
and n k ^ (0) = 0. 
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Figure 3: Velocity dependence of the diffusion constant D(V ) from the adiabatic to the 
diabatic limits. We use X — 1, fl/p — 0.1, poT^ = 2n (p thus defines the energy units), 
KqPo = 5 (the lowest curve) and K p = 15 (the highest curve). 
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Figure 4: Deviation from the fluctuation-dissipation theorem as a function of f3 — 1/T. 
For low velocities, the theorem is largely satisfied. However, for large velocities it is 
violated even at moderate temperatures. Here we use X — 1, p T^ = 2n and K p Q = 5. 

dependence (in the diabatic limit). At high velocities, the system becomes increasingly 
transparent, as reflected in the decrease of the diffusion constant. A similar behaviour 
is observed for the first cumulant, i.e. the average heating rate, which as a function of 
the velocity V has a similar profile with D(V ). This is reminiscent of the motional 
narrowing phenomenon in NMR. At low enough velocities the fast system has sufficient 
time to "accommodate" to the new environment, while the shape X changes. In the 
opposite limit of high velocities, the shape X evolves so rapidly that the system can 
barely react to the changes. Consequently, the energy diffusion is maximal only for some 
intermediate velocity regime, when the "slow" motion is in "resonance" with the "fast" 
dynamics, namely when t s1ow = X /V is comparable to Tf ast = Ti/kq. In Fig. 4, we plot 
the ratio of PD/j, where D and 7 are computed from the first and second order cumulants 
as a function of velocity V . When the ratio is unity, the Einstein limit of the fluctuation- 
dissipation theorem is recovered. Noticeable differences occur at large velocities. Higher 
order cumulants have also been extracted from our numerical results and in all cases their 
magnitudes and temporal behavior is similar to the analytical results discussed in the 
previous two subsections. Higher order cumulants increase essentially linearly with time 
and their magnitudes increase with the order of the cumulant. 
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4 Dynamical evolution of a simple quantum system 
in a complex environment 



In the adiabatic approximation (discussed in the previous section), it is possible to an- 
alytically determine the influence functional. Since the characteristic time scale for the 
intrinsic subsystem Tf ast = H/kq is much shorter than any expected time scale of the 
"slow" subsystem in Eq. (|T^) one obtains significant contributions only from "left " Y(t") 
and "right" X(t') paths corresponding to t" = t', and the influence functional acquires 
the particularly simple form given in Eq. (|i8|) . For a derivation of this form we refer the 
reader to Ref . |Tl| . As a result the evolution equation for the density matrix of the "slow" 
subsystem is local in time and memory effects are absent. Among the conditions we have 
listed for the applicability of our adiabatic results, we had j3 = 0, which corresponds to 
an infinite temperature of the intrinsic subsystem (since (3 = 1/T). As we have discussed 
in Section II, in the present formulation only heat transfer is allowed between the two 
subsystems. Since the temperature of the intrinsic system is infinite in this limit, one 
naturally expects energy transfer only from the intrinsic system towards the slow system 
and no mechanical work. Even though there are no memory effects in the time evolution 
of the "slow" subsystem, memory effects are still present in the evolution of the intrinsic 
subsystem (see Section III). 

It is worth noting that the functional form of the influence functional derived by us 
is entirely different from the popular Caldeira-Leggett form 0, which is a quadratic 
expression in X(t) and Y(t). If we were to use only the first term in a Taylor expansion 
of G(X(t'), Y(t')) — 1 we would obtain an expression similar to Caldeira-Leggett form for 
the influence functional. 

It is well known that most of the trajectories in a path integral are very jagged and 
one might question the applicability of an adiabatic approximation. One can try to 
make a long argument in support of the adiabatic approximation using "circumstantial" 
evidence, e.g. the success of the Born-Oppenheimer approximation in studying various 
physical systems. We only note that relatively smooth paths, obtained by taking into 
account only low frequencies, seem to be adequate, as is clearly shown in the partial 
Fourier smoothing methods for computing path integrals |18| . 

In this section we shall discuss several cases: that of a linear potential, a quadratic 
potential and a double well potential where dissipative tunneling can occur. Surprisingly 
enough, the first two cases can be solved explicitly. To our knowledge, these represent new 
cases of time dependent quantum mechanical problems, which can be given an entirely 
analytical treatment p3l. Only the case of a double well potential is treated numerically. 
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4.1 Linear potential 

The case when the slow variables evolve in a linear potential 

H " (x ^-&fi- FX 



(33) 



is particularly instructive. We shall assume furthermore that the "effective potential" is 
"translation invariant", namely G(X, Y) = G(X — Y). In terms of the variables s = X — Y 
and r = (X + Y)/2 the evolution equation for p(X, Y, t) becomes 



h 2 

(ihd t + — d r d s )p(r, s, t) = {-Fs + iT l [G(s) - l]}p(r, s, t). 
m 



We seek a solution in the form 

p(r,s,t) 



dk 
2~tx% 



exp 



/ ikr \ 

V*) 



d(k, s, t). 



(34) 



(35) 



The function d(s, t, k) satisfies the equation 



d t + — d s j d(k,s,t) 



— + -[G(s)-l]\d(k,s,t). 



(36) 



For either s = or k = 0, d(k, s, t) is the characteristic function [|T^] for the spatial or 
momentum distribution of the slow subsystem, respectively 

Using the method of characteristics for wave equations Eq. (36) can be solved 
through quadratures and the density matrix is determined to be 



p(r,s,t) 



exp 



x exp 



dr'dk ( . kt" 

— — p \r,s 

zirn \ m t 

'iFst iFt 2 k T l m 
+ 



ikir — r') 

n 



2hm 



hk 



(37) 



where po(r, s) = p(r, s, 0) is the initial density matrix. From Eqs. (j35[) and (]37f) , the char- 
acteristic function for the momentum distribution can be identified as D(s, t) = d(0, s, t), 
where: 



D(s,t) = J drp(r,s,t) = J drpo(r, s) exp 



(3? 



One extremely economical and intuitive way to characterize the momentum distribution 
of the collective subsystem is through its cumulants 



(39) 



s=0 
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For the case of a Gaussian correlation function G(X) = exp[— X 2 /2Xq], we find 



t=o + Ft, (40) 



((P 2 ))\t = <(P 2 »I*=0 + — , (41) 



VI ( h N 



((P 2n )>| t = (( P 2n ))\t=o + (2n - lyr— ^-j (42) 
((P 2n+1 ))k = ((P 2n+1 ))\t=0. (43) 

The meaning of the "correlation length" X is that intrinsic shapes separated by \X — Y \ > 
X are statistically uncorrelated. Notice that only the first cumulant is affected by the 
presence of a linear potential in the expected manner, namely a uniform acceleration of 
the slow subsystem. The "bath" of intrinsic degrees of freedom affects only higher order 
even cumulants of the momentum distribution while the odd cumulants of order higher 
than one remain unchanged. 

The cumulants of the spatial distribution can be obtained from the characteristic 
function d(k,0,t). From Eq. (|37|) it immediately follows that 

lnd o (A;,0,t) + 





+ -wL ds ' [G ^- l] - (44) 



The term In do(k, 0, t) gives the contributions to the cumulant expansion arising from free 
expansion of the initial wave packet, in the absence of both the linear potential and the 
coupling to the internal degrees of freedom. The linear potential leads to the expected 
(classical) behaviour of the center of the wave packet (see the second term on the rhs of 
Eq. (f44|) ). The contribution to the even cumulants arising from dissipation alone is 



ir 2n ))\ diss = { \^ V (z^l • (45) 
Of particular interest is the second cumulant 

"i- = S?< (46) 

which shows that dissipation leads to a super diffusive expansion of the wave packet. This 
behaviour is to be contrasted with the free expansion or ballistic propagation, in which 
case ((r 2 )) oc t 2 and with normal diffusion, for which ((r 2 )) oc t. It is instructive to estimate 
also the "size" of this state in phase space. In the limit t —>■ oo the dissipative contribution 
dominates, and one has 

ArAp w [((r 2 » \ dlss ((p 2 )) {#.,] = ^ 2 ■ (47) 
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Thus the "blob" spreads in phase space at a much faster rate than in the case of a simple 
quantum mechanical system. 

Similarly, one can show that with logarithmic accuracy for large values of the variable 
6{r,t) = \r — Ft 2 /2m\/t (assuming vanishing initial average linear momentum) the spatial 
distribution behaves as 



p(r, 0, t) oc exp \-u6(r, t) ln 1/2 6(r, t) 



(4£ 



where v is some constant. Thus the effect of dissipation is undeniably not only significant, 
but also leads to qualitatively new features. 



4.2 Quadratic potential 

Another case that is susceptible of an analytical treatment is that of a quadratic potential 
for the collective subsystem 

h 2 „ 9 muj 2 X 2 , , 

H (X) = -—d\ + ——. 49 
2m 2 



Using the representation for p(r, s, t) of Eq. fl35|), the equation for the transformed density 
now becomes 

(d t + —d s -muj 2 sdk] d(k,s,t) = ^r[G(s) - l]d(k,s,t). (50) 
\ m J n 

The method of characteristics |H| can again be used to determine its solution 



. f dk 

P(r,s,t) = y^-exp 



ikr 



x do \ s cos cut sin oot, mus sin oot + k cos cut ] (51) 

y muj J 

{pi /•* ~) 
— / dt'[G(s cos u(t-t') sinu(t-t')) - 1] }, 
n Jo mui J 

where do{s,k) = d(s,k,0). In a similar manner to the one described in the previous 
subsection, one can determine various cumulants. For both spatial and momentum dis- 
tributions only even cumulants are affected by dissipation |22| 

rl ( % \ 2n r wt 

((p 2n ))Uss = (2n-l)!!— f — J J o drcos 2n r (52) 
„ [(2n-l)!!F Vk (h\ 2n 



huj \ muXo / 
\{2n- l)!!] 2 Th ( h 
2 n n\ h \muX 
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2n 



(55) 



There is a noticeable difference with the case of a linear potential, in that all cumulants 
increase now only linearly with time. It looks as if the quadratic potential has a "focusing" 
effect on the spatial distribution. The fact that the momentum and spatial distributions 
are so similar should come as no surprise in the case of a harmonic oscillator, which 
possesses an obvious symmetry between the momenta and coordinates. In the limit t — > oo 
we obtain that the "uncertainty relation" has the following expression 



1/2 irhW t 
muiXQ 



ArAp « ((r 2 ))\ dlss ((p 2 ))Uss = (56) 



which has to be contrasted with the quadratic time behaviour in the linear or no potential 
cases, see Rel. (|48|). 

It is a simple matter to analytically continue these expression to the case of an inverted 
parabolic potential or barrier. This rather innocuous procedure, leads however to an 
entirely different time dependence of the cumulants, all of them increasing exponentially 
with time in this case (as cost and sinr become coshr and sinhr respectively). 



4.3 Tunneling in a symmetric double well potential 

The double well potential we analyse has the form 

For a strong enough barrier the spectrum of a relatively large number of low lying eigen- 
states is made up of doublets in the absence of the coupling to the intrinsic subsystem. 
The corresponding eigenfunctions are approximately symmetric and antisymmetric com- 
binations of wave functions localized in the two wells. The first nine eigenvalues are shown 
in Fig. 5 for the case a = 1/2 and 6 = 5 together with the potential. We have chosen as 
an initial state a linear combination of the eigenfunctions of the ground and first excited 
states (n = 0, 1) 

H (X) = -^-d 2 x + V(X), (58) 
H Q (X)MX) = E n ip n (X) (59) 

p(X,Y,0) = ^[MX)-MX)][MY)-MY)]. (60) 

In this case the particle is initially localized in one well with a probability almost equal to 
one, and only an exponentially small amount is present in the other well. In the absence 
of coupling to the intrinsic subsystem this state will tunnel almost entirely to the other 
well in an exponentially long time T tunne i = irh/(Ei — Eq), since the splitting between the 
two states is exponentially small. For the particular choice of parameters we have chosen 
(ft = l, m = 1/2, a = 1/2 and b = 5) r tunnel « 2, 300. 
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Figure 5: The first nine eigenvalues for the Schrodinger equation with the potential 
(|57|) with h — 1, m — 1/2, a = 1/2 and 6 = 5. (The splitting of the first two levels is 
approximately 1CT 3 .) 
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The coupling to the intrinsic subsystem is characterized by two parameters and Xo, 
which we have varied independently. We retain a Gaussian correlator G(x) = exp(— x 2 /2). 
As we show in Figs. 6 and 7, the effect of dissipation on the tunneling process is profound. 
Not only the tunneling rate is changed by orders of magnitude, but the shape of the wave 
packet is qualitatively different from the one in the absence of coupling. In the case of 
an usual quantum mechanical tunneling, the wave packet would have simply gradually 
changed from a state localized in one well to a state in the other well. Since the potential 
well is symmetric, the shape of the wave packet after a time T tun nei = it%/{Ei — E ) would 
have been simply the mirror image with respect to the origin of the initial one, see Fig. 



The role of the spreading width is rather simple to understand. Since the coupling 
between the two subsystems is defined by it, see Eq. (§), it is natural to expect that with 
increasing the rate of tunneling increases as well. We observe an approximate power 
law relationship between the probability to find the particle in the other well at a given 
time 



and the spreading width T^- (see Fig. 7), namely P+(t) oc T^ a , where a ~ 3/2. 

The role played by the correlation length Xo is however somewhat more subtle. The 
rate of tunneling is increasing dramatically when X is decreasing. At the formal level 
this can be easily understood from the evolution equation (|19|) for the density matrix. 
In the Xy-plane when \X — Y\ < Xo, the term responsible for dissipation is negligible. 
The smaller the correlation length Xo the bigger is the area where dissipation is directly 
effective. In the region \X - Y\> X the term iT l [G(X,Y) - 1] « -iT l dominates the 
time evolution of the density matrix. Thus, in the long time limit the density matrix 
tends to become almost diagonal. This tendency could have been inferred as well from 
the time dependence of the cumulants for the momentum distribution, which we have 
obtained in the previous two subsections. Since these cumulants increase with time, that 
means that the derivatives of the density matrix p(X, Y, t) with respect to s = X — Y 
become larger and larger for X = Y. 

Another and perhaps a better way to understand the role played by the correlation 
length X on tunneling is by returning to the initial picture, see Fig. 1. With decreasing 
X the intrinsic subsystem undergoes more and more transitions on its way from one well 
to another, since the density of (avoided) level crossings per unit length increases. We 
remind the reader, however, that trying to understand the role of the energy exchange 
between the two subsystems in terms of isolated jumps at level crossings is not quite 
correct, as it was discussed at length in Ref. ||. This also follows from the results 
presented in Section II. As we stressed there, had the picture of a random walk in the 
energy space been an appropriate one, the emergence of rather long tails in the energy 
distribution would have been rather difficult to explain. 

As our numerical results suggest and the structure of the evolution equation (|iD|) also 



6. 




(61) 
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Figure 6: The density profiles for t — 0, 1, 2, 3, 4 and 5 for the initial condition 
defined in Eq. (|60|) , in the double well potential of Fig. 5 for two different values of 
T^. The dashed line shows the density profile in the absence of dissipation after a time 
t = r tunne; = nh/iEh - E ) w 2, 300. 
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1.000 
0.500 




Figure 7: The probability to find the particle in the other well has the approximate 
power law behaviour P+(t) oc T^ 2 . Different curves correspond to times t — 1, 2, 3, 4 
and 5 respectively in ascending order. 
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seems to confirm, for t — > oo the density matrix has the limit p(X,Y,t) — > 5(X — Y), 
which is likely the only stationary solution of Eq. (|19|). At some point in time however, 
the adiabatic approximation we have used to derive the evolution equation (|T9| ) will 
cease to be valid, namely when the characteristic momenta will become of the order of 
mXQKo/h. In Section III we established that for velocities of the order of XqKq/K one has 
the crossover from the adiabatic to the diabatic regime, when the "motional narrowing" 
sets in. For higher velocities dissipation becomes weaker, though not at all negligible. 
At this point there is however one essential new element which enters the evolution of 
the "slow" subsystem, the memory effects. We have commented on this previously in 
the text. In the absence of the memory effects in the dynamics of the "slow" subsystem 
"thermalization" between the two subsystems apparently cannot be achieved. 



5 Conclusions 

We are still quite a distance from solving the problem we have set out to do in the 
introduction. Nevertheless we have hopefully succeeded in clarifying several aspects. Here 
is the best place to draw the line and evaluate what we have achieved so far and what 
remains to be done. 

We have developed a formalism, which describes a relatively simple quantum me- 
chanical system coupled to a "bath" of intrinsic excitations. We solved the dynamical 
evolution equations for these systems, without making any uncontrollable approximations 
or assumptions. We have analysed two types of problems, each of which being interesting 
in its own right and also relevant for understanding various types of experiments. The 
first class of systems are the so called driven systems, when a certain number of externally 
controlled parameters are changed. Systems ranging from complex molecules to quantum 
dots in variable magnetic or electric external fields can be studied in this way. One ques- 
tion which we hope we have shed light on, is the velocity dependence of the diffusion 
constant. During the years there have been quite a range of answers to this question, 
some of them rather intriguing . 



There are situations, when such parameters become dynamical variables, as for exam- 
ple in nuclear fission, for which one has to use the second type of methods described above. 
Our approach is based on an almost entirely microscopic description of the intrinsic sys- 
tem, using parametric random banded matrices. We have established that the dynamics 
show some new features. Perhaps the most prominent one is the appearance of extremely 
long tails in either energy, momentum or spatial distributions for the subsystems in in- 
teraction and the manifestly non-Gaussian character of the dissipative dynamics. These 
features raise significant doubts concerning the applicability of various phenomenologi- 
cal transport approaches, such as Fokker-Planck and Langevin equations ||25|| , to finite 
many-body quantum systems. 

Refinements of details of the present scheme are desirable, for example, the partic- 
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ular parametrization suggested in Ref. |13|] for the correlator (H) and a more realistic 
parametrization for the energy dependence of the average density of states p(e). The 
introduction of a correlator that is not translationally invariant, which is well inside the 
scope of our formalism, allows one to tackle the problem of the scattering by a localized 
complex system. 

We do not have yet a solution for the two subsystems in interaction outside the adia- 
batic limit and consequently, we do not have an understanding of the dynamical evolution 
for this very important for experiments case. The main difficulty here seems to reside in 
finding a practical scheme to solve the evolution equations, which are already known. An 
extension of the present formalism to open quantum systems (to include particle emis- 
sion) is desirable as well. Apparently, from the formal point of view, the only new element 
would be the replacement of the Hermitian Hamiltonian Hi(X) by a non-Hermitian one 
and that would not lead to unsurmountable difficulties. A non-Hermitian H\(X) would 
lead however also to a formalism in which the probability is not conserved anymore, but 
alternative approaches can be also envisioned. 

One limitation of our approach is the neglect of the shape dependence of the aver- 
age density of states for the intrinsic subsystem p(e,X). This limitation excludes the 
description of processes in which mechanical work is exchanged. There is no apparent 
technical or methodological difficulty in taking this additional feature into account and 
its implementation should be straightforward. Dissipative effects arising from the shape 
dependence of the average density of states alone, i.e. due to the time dependence in 
p(e,X(t)), are expected to give rise to the so called one-body type of friction[|(| or the 
wall-formula. A simple example is an ideal gas in a container, which changes its shape. 

We have not explicitly specified the number of degrees of freedom of the "slow" sub- 
system. The reader might be left with the impression that our approach is limited to 
one dimensional "slow" motion only. Actually, that is not the it would be clear 

from a closer analysis. There is however an element of the dynamics, which would appear 
only when the dimensionality is two or higher, namely the emergence of abelian and/or 
nonabelian gauge fields |L5| [U|]. By their definition, gauge fields couple to momenta and 
thus are expected to appear in the next to the leading order of an adiabatic approxima- 
tion. Naively one would incorrectly conclude that gauge fields are therefore non-adiabatic 
in nature. Gauge fields would lead undoubtly to new and interesting dynamical effects. 
The generalization of the present approach to include them seems to raise no technical or 
methodological issues. 

There always is the nagging question about a microscopic justification of the entire 
parametric random matrix approach and especially about trying to get an idea about 
where different key parameters come from, e.g. Xq,Kq,TK Obviously, that question 
lies outside the scope of this work. We think, however, that one can get a pretty good 
idea about the appropriate values for such parameters from both many-body theoretical 
models and experiment as well. For example, one can easily consider an ensemble of 
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many interacting Fermions in various continuously changing containers and thus extract 
Xq. We do not see a principial obstacle here. Concerning the applicability of the very 
idea of a random matrix, in order to model a many-body system, that seems to have been 
answered convincingly |T], |2|, [3], f|, [5| . 

A more ambitious goal is to extend the formalism as to attain a more complete char- 
acterization of the two subsystems in interaction. So far the formalism conserves only 
the probability. One might want to have a formalism in which the total energy or other 
relevant observables, such as linear and/or angular momentum, is conserved. For exam- 
ple, one can devise evolution equations for a slightly generalized density matrix, such as 
p(X, Y, e, t), which is obtained by summing over intrinsic states with a given energy e only. 
This aspect has been barely touched upon by us p7] and appears to be a most daunting 
task. One can write down with relative ease an impressive number of evolution equations, 
satisfying such requirements. However the analysis and the solution of these equations 
seems at the present moment quite difficult, to say the least. Another direction, which 
seems to be the most difficult technically at the present time, is to extend the formalism 
so as to be able to compute fluctuation characteristics of the dynamics. 

We greatly appreciate the DOE support for this work and the computing facilities pro- 
vided by IDRIS and NERSC. The Laboratoire de Physique Theorique et Hautes Energies 
is a Laboratoire associe au C.N.R.S., URA 0063. 
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